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PREFACE 


This  Memorandum  discusses  the  motion  of  electrons  in  a  particu¬ 
lar  radiofrequency  field  to  determine  whether  confinement  is  likely 
in  this  environment.  The  problem  is  studied  for  two  cases:  one  in 
which  collisions  of  electrons  with  other  particles  are  negligible  (the 
collisionless  case),  and  one  in  which  electrons  collide  only  with 
neutrals  (the  weakly-ionized  case).  In  the  first  instance,  the  motion 
of  the  electrons  is  such  that  they  will  always  remain  in  the  interior 
of  a  certain  region  determined  by  the  value  of  the  magnetic  field, 
initial  electron  velocity,  and  tube  radius,  as  long  as  xhe  nagaetic 
field  strength  and  field  frequency  satisfy  certain  inequalities.  In 
the  second  case,  the  electrons  will  remain  within  the  tube  for  all 
values  of  magnetic  field  strength;  however,  since  diffusion  is  neglect¬ 
ed  in  the  basic  equations,  this  result  does  not  imply  that  the  elec¬ 
trons  are  confined. 

The  results  of  this  study  are  applicable  to  the  fields  of  electric 
propulsion,  controlled  thermonuclear  fusion,  and  gaseous  electronics. 
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SUMMAHY 


The  motion  of  electrons  in  a  radiofrequency  field  is  investigated 
for  a  cylindrical  geometry  and  an  electromagnetic- field  configuration 
consisting  of  a  spatially  constant  axial  magnetic  field  and  an  azi¬ 
muthal  electric  field  which  varies  linearly  with  distance  from  the 

r 

center  of  the  cylinder.  This  field  configuration,  which  approximates 
that  found  in  the  electrodeless  ring  discharge  near  breakdown,  has 
also  been  suggested  as  a  plasma -confinement  scheme.  When  the  modioli 
is  collisionless,  the  equations  can  be  solved  analytically  to  show 
that  stable  (time -bounded)  motion  exists  for  certain  values  of  the 
ratio  of  electron  cyclotron  frequency  to  applied  field  frequency.  Con¬ 
finement  of  the  orbits  to  a  region  entirely  inside  the  cylinder  is  com¬ 
pletely  assured  only  for  those  particles  which  start  from  rest;  in  all 
other  cases  the  maximum  radius  of  extent  depends  on  initial  position 
and  velocity. 

Collisions  are  quali tat i vely^ac counted  for  in  a  crude  manner  by 
use  of  the  Lange vin  equation,  which  includes  a  drag  force  due  to  colli¬ 
sions  in  the  equations  of  motion. ftvEven  with  this  simple  model  the  re¬ 
sulting  equation  is  apparently  intractable  except  by  machine  solution. 
It  is  therefore  not  possible,  without  analytic  solutions,  to  discuss 
the  effect  of  collisions  of  this  type  on  the  stability  of  the  motion, 
although  it  appears  from  several  of  the  examples  in  this  Memorandum 
that  an  initially  unstable  or  divergent  motion  quickly  damps  out  when 
the  friction  force  is  added,  and  that  for  large  friction  forces,  the 
mean  motion  is  always  stable.  Confinement  is  not  ensured,  however, 
because  diffusion  is  not  included  in  the  basic  equations. 

Some  numerical  examples  of  the  motion  are  presented  for  both  the 
collisionless  case  and  the  Lange  vin  equation.  An  approximate  solu¬ 
tion  is  obtained  for  the  case  with  collisions  which  adequately  de¬ 


scribes  the  mean  motion  whenever  the  ratio  of  collision  frequency  to 
field  frequency  is  larger  than  three. 
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v  =  auxiliary  function 

X  =  ~'(Eq.  (B-l) ) 
y  =  2  Re( p)  *  y  sin  t  (Eq,  (B-8)) 

z  =  position  coordinate 

a .  .  =  coefficients  in  general  solution  (Eq.  (A-6)) 

^  J 

0  =  ratio  of  collision  frequency  to  field  frequency 

y  =  Q/p,  ratio  of  cyclotron  frequency  to  collision  frequency 

0  =  angular  coor dilate 

\  =  real  parameter  (Eq.  (A-2)) 

=  K  +  ia  =  k/2 n,  characteristic  exponent 
v  =  colli son  frequency 

p  =  (v  Sin  T  +  i)/2  (Eq.  (B-l)) 
t  =  u>t,  dimensionless  time 

$  =  periodic  function  of  T  (Eq.  (A- 11)) 

p 

cp  =  Cl  sin  t/2 

n  =  ratio  of  cyclotron  frequency  to  field  frequency 
a)  =  field  frequency 
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I,  I  IfTRODUCT  IOH 
( 1-7 ) 

Several  recent  papers  have  been  concerned  with  the  confine¬ 

ment  of  plasmas  by  radiofrequency  fields.  Some  of  these  treat  the 
problem  from  the  continuum  standpoint;  others  have  used  the  single - 
particle-orbit  theory.  In  the  latter  category,  Weibel^  has  investi¬ 
gated  the  radial  bounds  for  electrons  in  a  TE^  mode  waveguide  at 
cutoff  and  has  concluded  that  the  orbits  are  stably  bound  for  certain 
values  of  the  ratio  of  cyclotron  frequency  to  applied  field  frequency. 

The  equations  of  motion  in  Weibel's  study  were  complicated  by  the 
electromagnetic- force  terms,  which  are  Bessel  functions  for  this  wave¬ 
guide.  He  investigated  the  solution  in  two  regions  where  the  equations 
became  tractable:  one  where  the  field  frequency  is  much  larger  than 
the  cyclotron  frequency,  and  another  for  small  radii,  valid  in  regions 
close  to  the  waveguide  axis.  In  neither  case  was  the  actual  motion 
found,  nor  were  the  effects  of  initial  conditions  discussed.  Further, 
in  the  second  case  an  apparently  unnecessary  assumption  was  made  in 
the  analysis,  which  yielded  an  erroneous  angular -momentum  equation, 
(ihis  assumption  did  not  affect  the  radial  motion. ) 

The  present  analysis  will  consider  the  motion  of  electrons  for 
a  situation  similar  to  Weibel's  second  approximation.  A  more  complete 
investigation  will  be  made  of  the  trajectories,  the  effects  of  boundary 
conditions,  and,  qualitatively,  the  effects  of  collisions  on  the  radial 
bounds  of  the  motion. 
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II.  THE  EQUATIONS  OF  M3TION 

The  electromagnetic  field  for  the  TE^  waveguide  of  Ref.  6  is 
given  by 


E9  ■  '  BoJl(f)CO,rt 

B  »  B  J  (  — ^  )  sin  cut 
z  o  o  \  c  / 


a) 


where  r,  9 ,  z  are  cylindrical  coordinates,  w  is  the  field  frequency, 
and  c  is  the  velocity  of  light.  The  Weibel  approximation  (cur/c)2  «  1 
is  equivalent  to  replacing  the  Bessel  functions  by  their  values  for 
small  argument 


E 


9 


r  cos  cut 


B  »  B  sin  cut 
z  o 


(2) 


One  can  visualize  this  approximation  as  a  restriction  to  small 
radial  distances,  which  is  necessary  if  the  frequencies  are  in  the 
microwave  regime.  On  the  other  hand,  if  the  frequencies  are  much 
lower  than  microwave,  the  radius  need  not  be  restrictively  small* 

This  particular  interpretation  of  Welbel's  approximation  permits  the 
coupling  of  a  plasma-confinement  scheme  to  an  operable  discharge  de¬ 
vice  of  current  interest — the  electrodeless  ring  discharge. 

A  cross-sectional  view  of  an  electrodeless  discharge  is  sketched 
in  Fig.  1,  in  which  an  axial  magnetic  field  is  induced  by  an  excita¬ 
tion  coil  of  several  turns,  and  a  Faraday  cage  can  be  placed  around 
the  tube  to  Isolate  the  interior  from  the  axial  electric  field  of  the 
coll.  When  the  discharge  is  operating,  the  currents  in  the  gas  give 
rise  to  cm  Intense  circumferential  glow,  which  accounts  for  the  "ring" 
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Fig.  I— Geometry  of  the  electrodeless  ring  discharge 
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terminology;  the  use  of  the  word  "electrodeless"  differentiates  this 
case  from  that  discharge  where  the  electrodes  come  into  contact  with 
the  gas. 

At  breakdown,  just  as  the  discharge  becomes  self-sustaining,  the 

(8) 

electromagnetic  fields  in  the  interior  are  given  by  Eq.  (2),  '  as 

the  electron  density  increases,  the  induced  currents  in  the  interior 

(a) 

distort  these  fields,  and  their  determination  can  be  quite  difficult. 
Therefore,  the  use  of  Eq.  (2)  in  this  analysis  precludes  large  electron 
concentrations,  and  the  results  will  apply  to  a  gas  in  which  the  elec¬ 
tric  conductivity  is  small.  A  treatment  of  the  case  in  which  the 

(7) 

densities  are  higher  is  given  by  Zauderer. 

To  obtain  the  motion  in  the  case  of  a  weakly  ionized  plasma,  we 
will  use  the  Lange vin  equation.  This  equation  has  been  used  extensi¬ 
vely  to  describe  the  ac  conductivity  of  a  plasma  (see,  e.g.,  Ref.  10) 
when  the  only  important  interactions  are  those  occurring  between  elec¬ 
trons  and  neutrals. 

From  Allis,  the  Lange  vin  equation  is  written 

U  -  ?  +  A(t)  (3) 

where  F  is  the  Lorentz  force  on  the  electron 

cw 

and  A(t)  is  a  fluctuating  force  due  to  the  collisions.  In  order  to 

solve  Eq.  (3),  the  force  A(t)  must  be  defined.  If  we  assume  that 
— »  — •  — * 

A(t)  is  independent  of  v,  and  that  A(t)  varies  extremely  rapidly  com¬ 
pared  to  the  variations  in  v,  then  it  is  possible  to  define  A(t)  sta¬ 
tistically  and  solve  Eq.  (3). 

— * 

The  second  assumption  made  for  A(t)  is  tantamount  to  assuming 
that  the  motion  of  the  electron  is  taking  place  in  two  ways;  that  is, 
there  exists  a  small  time  interval  At  during  which  the  variations  in 
v(t)  are  small  while  A(t)  may  oscillate  several  times.  Furthermore, 
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since  l(t)  is  a  random  force,  there  is  no  correlation  between  A(t) 
and  A(t  +  At).  (Such  a  process  is  called  a  Markoff  process, in 
which  the  occurrence  of  an  event  at  time  t  depends  only  on  the  state 
of  the  system  at  time  t. )  Even  so,  the  velocity  defined  by  Eq.  (3) 
must  vary  in  a  continuous  way  so  that  a  solution  of  Eq.  (3)  will  have 
meaning  in  a  physical  sense.  This  smoothing  process  is  a  serious 
assumption,  and  the  use  of  the  Lange vin  equation  in  this  particular 
case  has  been  suggested  only  because  of  the  a  posteriori  success  of 
this  equation  in  predicting  with  reasonable  accuracy  certain  trans¬ 
port  phenomena  in  electromagnetic  fields.^11*’ 

Because  the  Langevin  equation  yields  the  motion  of  an  average 
electron,  i.e.,  a  drift  velocity,  it  is  more  suitable  for  the  determi¬ 
nation  of  mobility  and  conductivity  than  for  the  present  study  of  con¬ 
finement,  where  we  are  concerned  with  the  loss  of  individual  electrons 
from  the  plasma.  The  inability  of  the  Langevin  to  produce  this  in¬ 
formation  occurs  because  diffusion,  which  is  a  primary  mechanism  for 
loss  in  any  plasma,  cannot  be  taken  into  account  in  this  "smoothed" 
picture  of  the  motion.  What  the  Langevin  equation  will  show,  diffu¬ 
sion  excluded,  is  how  the  average  electron  moves  about  in  the  con¬ 
tainer  under  the  influences  of  the  electromagnetic  forces  and  colli¬ 
sions  with  neutral  particles.  Naturally,  if  this  motion  is  toward 
the  container  walls,  or  if  it  grows  large  with  time,  one  would  suspect 
that  the  forces  are  not  favorable  for  confinement.  However,  because 
diffusion,  reattachment,  and  recombination  are  neglected,  one  cannot 
infer  that  a  stable  Langevin  orbit  implies  confinement. 

With  the  understanding  that  ^(t )  must  have  a  large  frequency 
compared  to  the  variations  in  v(t),  then,  we  can  define  the  fluctua¬ 
tion  term  as 


A(t)  *  -  v v  (5) 

where  v  is  a  collision  frequency  which  must  by  assumption  be  higher 
than  any  of  the  other  frequencies  in  the  equations  of  motion.  It  can 
be  identified  with  the  collision  frequency  for  momentum  transfer  vc 
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if  the  particle-distribution  function  is  Maxwellian.  In  this  case  v 
is  independent  of  the  energy  of  the  particle,  and  the  first  assumption 
for  A( t)  is  met.  Molmud^1^)  has  investigated  the  solution  of  Eq.  (3) 
when  F  *  EQeiu/t  and  has  shown  that  when  the  distribution  function  is 
non-Maxwellian,  v  may  be  complex  and  a  function  of  temperature.  When 
I v |  «  a),  or  j v |  »  u),  the  imaginary  part  of  v  is  negligible.  We  will 
restrict  ourselves  to  real  values  of  v  in  this  analysis. 

Combining  Eqs.  (2)  through  (5)  gives  the  basic  equation  of  motion 
for  the  system  under  consideration: 


The  forces  in  Eq.  (2)  are  such  that  the  motion  of  the  electron 
is  entirely  in  the  cross-sectional  or  rS  plane.  It  is  convenient  to 
write  Eq.  (6)  in  complex  form,  letting  z  *  re10.  One  obtains 

zf  f  -  i  ^  T  2  +  (Q  sin  t  +  i0)z’  J  ■  0  (7) 


where  z  has  been  made  dimensionless  with  respect  to  the  initial  posi¬ 
tion  r(0),  the  primes  denote  differentiation  with  respect  to  time 
t  *  art,  and  the  dimensionless  parameters  0  and  0  are  defined  as 


eB 

_ o  cyclotron  frequency 

mcu;  field  frequency 


and 


v  m  collision  frequency 
0)  field  frequency 


(8) 


(9) 


The  initial  conditions  in  Eq.  (7)  are  z(o)  *  1  and  z'(o)  «  V  +  iC^, 
where  V  and  are  dimensionless  velocities  in  the  radial  and  angular 
directions. 
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Under  the  conditions  of  averaging  discussed  previously,  z  is  the 
position  of  an  average  electron  at  time  t  when  3  is  larger  than  the 
associated  frequency  of  z',  Before  discussing  the  formal  solution  of 
Eq.  (7),  let  us  determine  the  restrictions  placed  on  3  and  Q  by  the 
Langevin  approximation.  Allis^11^  discusses  the  case  in  which  Eq.  (6) 
is  solved  for  a  constant  magnetic  field  and  a  sinusoidal  electric 
field.  Although  Allis  did  not  specify  the  value  of  the  time  average 
in  comparison  to  the  various  frequencies,  it  can  be  deduced  from  his 
treatment  that 


1  < 
3  ~ 


^  «  mAt  «  1 


(10) 


i.e.,  he  averaged  over  a  time  At  greater  than  the  cyclotron  period  but 
less  than  the  period  of  the  field.  Since  the  energy  is  constant  over 
a  cyclotron  period  in  a  static  magnetic  field,  the  energy  assumption 
is  satisfied,  and  the  velocity  varies  as  the  field  frequency  u>. 

On  the  other  hand,  Molmud^1^  discusses  the  same  situation  and 
treats  both  the  case  of  3  »  1  and  that  3  «  1.  On  the  basis  of  the 
discussion  by  both  Chandrasekhar^ 12  ^  and  Allis^11^  that  the  time  At 
must  be  long  compared  to  l/v  and,  furthermore,  that  the  energy  must 
be  constant  during  At,  it  is  difficult  to  see  the  significance  of  the 
case  of  3  «  1,  which  is  equivalent  to  the  inequality 


1 

uu 


« 


1 

V 


(11) 


or,  since  the  collision  time  must  be  shorter  than  the  averaging  time, 

-  «  At  (12) 

U) 

Ibis  inequality  indicates  that  the  average  has  been  taken  over  the 
period  of  the  applied  field,  or,  equivalently,  that  the  energy  is 
constant  over  the  period  of  the  applied  field.  This  averaging  process 
may  be  used  to  define  an  effective  electric  field^11^  which  is  inde- 
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pendent  of  the  periodic  terras  and  may  be  treated  as  a  dc  field  in  the 
breakdown  theory  for  an  ac  discharge.  In  the  present  case,  however, 
an  averaging  process  such  as  the  above  would  not  yield  meaningful  re¬ 
sults. 

Therefore,  solutions  of  Eq.  (7)  will  be  discussed  for  the  case 
given  by  Eq.  (10)  and  the  case  in  which  the  average  time  is  shorter 
than  the  cyclotron  period: 

|  «  3  ~  1  (13) 

We  will  also  obtain  solutions  of  Eq.  (7)  for  3=0. 

The  general  solution  of  Eq.  (7)  is  discussed  completely  in  Appen¬ 
dix  A.  For  convenience,  the  solution  and  subsidiary  equations  are  re¬ 
peated  below.  After  several  transformations,  one  obtains  for  the 
position  of  the  average  electron 

z  =  v(T)e{if'  sin2  7/2  -  3t/2)  (14) 

where  w  is  the  general  solution  of  the  differential  equation 

v  »  .  0  (15) 


with  initial  conditions 


w(0)  =  1 

w '  (0)  =  V  +  iC  +  3/2 


(16) 


Equation  (15)  reduces  to  the  Mathieu  equation  when  3*0,  vhich 
is  consistent  with  Weibel's  solution.  ^  When  3  ^  0,  a  general  solu- 
tion  is  more  difficult  to  obtain.  The  collisionless  motion  will  be 
discussed  first,  so  that  the  general  behavior  of  the  undamped  motion 
may  be  ascertained. 
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THE  COLLISIONLESS  MOTION 

Although  the  collisionless  case  has  been  previously  derived  by 
Weibel, ^  neither  explicit  trajectory  information  nor  initial  con¬ 
ditions  were  discussed.  Since  we  are  interested  here  in  both  the 
motion  and  the  possibility  of  containment  (compared  to  boundedness), 
we  will  investigate  Eqs.  (14)  through  (l6)  more  thoroughly  in  this 
section. 

When  0*0,  Eq.  (15)  becomes  the  Mathieu  equation,  which  is  dis¬ 
cussed  more  thoroughly  in  Appendix  A.  Let  X(t)  and  Y(t)  be  solutions 
to  Eq.  (15);  then  Eq.  (l4)  becomes 

R2  .  |z|2  .  (X  +  VY)2  +  Y2  (17) 


and 


0 


arg  z  a  tan 


ciY 

X  +  VY 


+  ft  sin2r/2 


(18) 


in  polar  coordinates,  where  R  is  dimensionless  with  respect  to  the 
initial  radial  position,  and  0(o)  »  0.  Because  the  collision  term 
is  zero,  it  is  not  necessary  to  average  Eq.  (3),  so  that  R  and  0 
describe  the  actual  position  of  the  electron. 

It  is  of  interest  to  compare  the  above  solution  with  the  analy¬ 
sis  made  by  Weibel.  If  one  solves  the  differential  equations  given 
in  Ref.  6,  which  are  written  in  Cartesian  coordinates,  it  is  seen 
that  the  radial  motion  is  identical  to  that  of  Eq.  (17)*  although 
Weibel  did  not  specify  the  initial  conditions.  On  the  other  hand, 
the  angular  motion  resulting  from  Weibel* s  analysis  lacks  the  periodic 
term  ft  sin2T/2,  due  to  an  assumption  made  by  Weibel  while  linearizing 
the  equations  of  motion.  Because  the  centrifugal-force  term  becomes 
singular  at  the  origin,  Weibel  saw  that  the  problem  would  become 
tractable  if  the  electromagnetic  forces  were  approximated  by  a  scalar 
potential,  which  removed  the  singularity  and  permitted  linearization 
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of  the  equations.  This  scalar  potential  affects  the  angular  motion, 
but  it  does  yield  values  of  the  radial  motion  which  are  correct.  Since 
Veibel  was  concerned  primarily  with  the  radial  bounds  of  the  orbits, 
the  conclusions  of  Ref.  6  are  not  affected  by  the  error  in  the  angular 
motion. 

The  singularity  still  exists  in  the  present  analysis,  since  both 
R*  and  0*  become  infinite  as  R  0  except  when  the  initial  angular 
velocity  CL  =0.  However,  this  singularity  will  never  be  attained; 
inspection  of  the  Hamiltonian  of  the  system^  ;  shows  that  the  conju¬ 
gate  angular  momentum  is  given  by 


po  5 


mr 


e*  - 


-)  sin 


(19) 


where  the  first  term  is  the  angular  momentum,  and  the  last  term  is  the 
azimuthal  component  of  the  vector  potential  for  the  field  configura¬ 
tion  of  Eq.  (l).  Because  we  defined  to  be  the  value  of  Q 1  when 
t  =*  0,  and  since  the  second  terra  vanishes  when  t  =0,  it  is  immediately 

obvious  that  when  C  =  0,  the  quantity  Pn  »  0  for  all  times,  and  the 

-L  y 

singularity  does  not  exist.  Similarly,  when  ^  0,  then  /  0  and 
the  orbit  will  never  pass  through  the  singular  point  at  R  ■  0.  This 
can  be  seen  in  Figs.  2  through  5» 

Therefore,  by  linearizing  Eq.  (l)  before  insertion  into  the  equa¬ 
tions  of  motion,  and  by  utilizing  the  initial  conditions  of  the  prob¬ 
lem,  it  is  possible  to  obtain  solutions  of  the  motion  even  when  the 
singularity  at  the  origin  is  not  removed.  This  refinement  is  not  im¬ 
mediately  obvious  if  one  works  in  cylindrical  coordinates  (the  natur¬ 
al  coordinate  system),  because  then  the  equations  are  nonlinear. 

The  behavior  of  R  and  0  is  obtained  by  investigating  the  proper¬ 
ties  of  the  Mathieu  functions.  The  Mathieu  equation,  as  discussed  in 
Appendix  A,  is  a  special  ca3e  of  the  Hill  differential  equation.  As 
such,  there  exist  regions  of  stable  and  unstable  solutions,  depending 
on  the  parameter  O.  The  regions  of  stability  are  tabulated  in  Ref.  15 
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and  plotted  in  Fig.  9  (Appendix  A),  where  it  is  shown  that  time-bounded 
solutions  exist  for  Q  in  the  following  regions:  0  <  Q  <  2.29;  3*78  < 
n  <  5.1+5;  6.96  <  0  <  8.60;  etc.  Thus  for  these  values  of  Q,  the  mo¬ 
tion  defined  by  Eqs.  (17)  and  (l8)  will  be  bounded  for  all  values  of 
time,  i.e., 


M  +  VY)2  +  C^Y2  <  M  (20) 

where  the  bound  M  is  sow  function  of  0,  V,  and  C^.  Itow,  in  order  for 
the  condition  given  by  Eq.  (20)  to  ensure  containment  of  the  particles 
within  a  given  region,  it  must  be  specified  further  that 

i<m<ro  (21) 

where  Rq  is  the  dimensionless  radius  of  the  tube. 

Unfortunately,  it  is  not  possible  to  determine  the  maximum  values 
of  |z|  for  a  given  set  of  initial  conditions  without  resorting  to 
numerical  calculations,  because  the  Mathieu  functions  X  and  Y  do  not, 
as  a  rule,  have  extensive  recursion  relations  between  themselves  and 
their  derivatives.  The  difficulty  of  obtaining  a  least  upper  bound 
for  |z|  is  illustrated  if  we  write  R  in  terms  of  the  actual  Jfethieu 
functions. 

For  0  in  the  range  of  stable  solutions  given  above.  Appendix  A 
shows  that  it  is  possible  to  have,  by  appropriate  choice  of  0,  peri¬ 
odic  solutions  of  Eq.  (15).  For  convenience,  let  us  choose  solutions 
so  that  the  Mathieu  functions  are  of  fractional  order;  when  the  frac¬ 
tion  is  rational,  say  a  3  p/s  <  1,  then  the  period  is  2rrs.  Hie  prob¬ 
lem  of  determining  a  for  a  given  0  is  discussed  at  length  in  Ref.  15 . 
Hie  Mathieu  functions  corresponding  to  these  parameters  are 
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X  3  cek(T,  Q) 


^  cos  (m  +  <x)-r 

ra=-»  m 


00 


Y  h  sek(T,  n) 


I  a! 


^  tS^  sin  (m  +  c*)t 
m=-»  in 


(22) 


1  4k) 

. _ _  m 


The  subscript  k  «  m  +  a,  where  mis  either  even  or  odd,  depending  on 
the  value  of  0. 

Inserting  Eq.  (22)  into  the  position  equations 


1,2  ■  K  *  v“k?  *  °?  4 


(23) 


and 


6  “  ^  ((ce^+  Vsek))  +  0  ein2j/2  (2U) 

From  Eq.  (2l)  we  can  write  the  restriction  on  R  as 

R2  -  [cek  +  VseJ2  +  C2  se2  *  l  s  M2  <  R^  (25) 

It  is  seen  from  Eq.  (25)  that  the  bound  M  on  R  will  be  a  function  of 
ft  (through  the  Mathieu  functions),  V,  and  Cy  The  straightforward 
procedure  of  differentiating  Eq.  (25)  to  find  the  of  R  in  the 

region  0  £  t  *  2ns  does  not  permit  a  solution,  because  of  the  lack  of 
analytic  relationships  between  ce,  se,  and  their  derivatives. 

Since  the  mathematical  determination  of  the  radial  bounds  is 
difficult,  let  us  examine  the  physical  situation  at  t  ■  0,  Starting 
from  R  ■  1  at  r  *  0,  the  only  force  acting  on  an  electron  is  the 
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azimuthal  electric  field,  which  will  cause  an  angular  acceleration 
which  decreases  in  magnitude  as  t  increases.  Suppose  the  electron 
had  an  outward  velocity  at  r  *  0  which  was  larger  than  the  azimuthal 
velocity  obtained  from  the  electric  field.  Then,  if  the  particle  is 
initially  near  the  outer  radius  (i.e.,  if  ^  0(l)),  the  trajectory 
could  intersect  the  wall  of  the  tube.  This  would  imply  that  the  only 
particles  always  contained  in  the  plasma  for  a  stable  solution  of  Eq. 
(15)  are  those  which  start  from  rest.  This  is  unlikely  in  a  gas  with 
nonzero  thermal  energy,  so  it  appears  that  some  particles  will  always 
be  lost  to  the  wall. 

Some  examples  have  been  prepared  to  illustrate  this  collisionless 
motion — particularly  the  effect  of  the  initial  conditions  and  cyclo¬ 
tron-frequency  parameter  Q  on  the  behavior  of  the  orbits.  The  solu¬ 
tions  were  obtained  on  an  IBM  yoyo  computer,  since  the  Mathieu  func¬ 
tions  desired  for  the  solutions  are  not  generally  available  in  tabu¬ 
lated  form.  (Earlier  solutions^ ^ ^  using  the  series  were  later  found 
to  have  convergence  problems:  in  particular,  the  simple  solution  for 
the  angular  motion  was  not  recognized,  and  the  numerical  integration 
for  the  angle  6  incurred  cumulative  errors.  In  some  cases,  however, 
the  series  solution  agreed  well  with  the  machine  solution,  as  will  be 
seen  in  Fig.  k . ) 

Figures  2  and  3  show  the  motion  for  a  natural  period  of  4n  in  the 
case  of  low  and  high  cyclotron  frequency.  Only  the  motion  for  2rr  is 
presented  in  Fig.  3;  during  the  second  half  of  the  period  the  electron 
retraces  its  path  back  to  the  starting  position.  There  is  one  value 
of  in  each  stable  region  of  the  Mathieu  equation  which  will  yield  a 
period  of  4n.  The  lowest  value,  0  a  I.387,  is  given  in  Fig.  2.  For 
a  discharge  at  30  me,  this  corresponds  to  a  weak  magnetic  field  on  the 
order  of  10  gauss. 

At  the  other  extreme.  Fig.  3  gives  the  motion  for  0  *  17.242, 
where  Bq  ~  175  gauss  at  30  me.  The  effect  of  increased  magnetic  field 
is  directly  evident  here,  where  the  electron  makes  5  circuits  in  a 
quarter-period.  The  order  of  occurrence  of  the  circuits  is  marked  by 
number  in  Fig.  3.  Notice  also  that  for  a  large  portion  of  the  trajec- 


Fig.  2  — Collisionless  motion  for  small  magnetic  field 


0  s  V  =  C,  =  0 
Period  Air 

A  Position  of  vorious  times 


Fig.  3  —  Collisionless  motion  for  large  magnetic  field 
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tory,  the  orbit  is  nearly  circular.  Here,  the  guiding- center  approxi¬ 
mation  is  valid, for  the  electron  can  make  several  gyrations 
in  an  essentially  constant  magnetic  field.  For  values  of  0  larger  than 
20,  the  guiding-center  approximation  could  be  used  for  that  portion 
of  the  trajectory  where  B  0. 

Figure  4  demonstrates  that  particles  with  outward  initial  velocity 
can  move  toward  the  container  wall  for  a  portion  of  their  trajectory. 
The  dimensionless  velocities  of  Fig.  1+  compare  to  a  thermal  velocity 
of  about  8V  (or  8C^)  ev  in  a  30-mc  discharge. 

Rather  complicated  trajectories  can  be  observed  when  a  particle 
with  initial  velocity  is  subjected  to  a  high  magnetic  field.  Figure  5 
shows  the  complex  motion  undergone  by  an  electron  during  one-half 
period  when  0  =  17.242  and  V  *  »  1.  Again,  we  have  indicated  the 

order  in  which  the  electron  traverses  the  tube.  For  a  stable  solution 
the  energy  gained  in  the  first  quadrant  must  be  lost  by  deceleration 
during  the  second  quarter  of  its  trajectory.  As  a  result  the  electron 
spends  a  considerable  portion  of  its  trajectory  in  a  small  region  of 
the  tube.  This  behavior  is  favorable  for  either  confinement  or  break¬ 
down.  During  the  reminder  of  the  period  the  motion  takes  place  pri¬ 
marily  in  the  third  and  fourth  quadrants. 

Although  this  discussion  of  collisionless  motion  has  been  brief, 
the  pattern  of  electron  motion  in  this  type  of  a  radiofrequency  field 
has  been  established:  Stable  solutions  exist  when  the  ratio  of  cyclo¬ 
tron  frequency  to  field  frequency  is  within  certain  limits.  When  0 
»  1,  the  motion  over  a  portion  of  the  trajectory  is  adequately  de¬ 
scribed  by  the  guiding-center  method.  If  ft  ±G  outside  the 

region  for  stable  solutions,  the  motion  is  divergent. 

An  example  of  divergent  motion  is  given  in  Fig.  6,  where  a  solu- 

(17) 

tion  originally  obtained  by  Brault  has  been  extended  for  larger 
times  to  indicate  the  increase  in  maximum  radius  as  t  increases.  The 
equation  derived  by  Brault  was  not  recognized  as  the  Mathieu  equation 
(in  cylindrical  coordinates  this  is  not  easy  to  see^^),  and  as  a 
result,  some  of  his  solutions  were  inadvertently  in  the  unstable  re¬ 
gion.  One  of  these  unstable  solutions  occurs,  as  evident  in  Fig.  9, 


Fig.  4 —  Collisionless  motion  with  initial  velocity 
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Fig.  5 —  Collisionless  motion  in  large  magnetic  field  with  initial  velocity 
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Fig.  6 —  Unstable  collisionless  motion 
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vhen  0  *  2.JL6;  this  is  the  case  plotted  in  Fig.  6  for  which  the  ini¬ 
tial  conditions  were  taken  as  V  ■  *  0.0205*  It  can  be  seen  that 

the  radius  of  gyration  increases  every  time  the  electron  passes  near 
the  origin  (because  V  *  ^  0,  the  trajectory  never  passes  through 

the  origin,  but  since  V  «  «  1,  the  distance  of  the  path  from  the 

origin  is  small).  The  kinetic  energy  also  exhibits  an  undamped  in¬ 
crease,  because  each  time  the  particle  passes  near  the  origin  it  gains 
energy,  as  compared  to  the  stable  motion  shown  in  Figs.  2  through  5, 
where  a  periodic  gain  and  loss  of  energy  is  achieved  as  the  electron 
traverses  the  tube. 

MOTION  WITH  COLLISIONS 

When  9^0,  Eq.  (ik)  gives  the  position  of  an  average  electron 

in  the  tube,  rather  than  the  motion  of  individual  electrons  as  it  did 

in  the  preceding  section.  This  occurs  because  the  collision  term  in 

Eq.  (6)  has  been  averaged,  and  while  momentum  is  conserved  during 

collisions,  no  changes  in  direction  are  made  as  a  result  of  collisions; 

hence  if  one  were  to  follow  an  individual  electron  as  we  did  in  the 

last  section,  the  position  given  by  Eq.  (l4)  would  not  describe  the 

correct  motion.  If,  however,  one  were  to  measure  a  current  j  *  nev  . 

ave 

Eq.  (14)  would  indicate  the  location  and  behavior  of  this  current 
resulting  from  the  group  behavior  of  electrons  in  this  radiofrequency 
field. 

Keeping  this  in  mind,  let  us  discuss  briefly  the  solution  of  the 
system  of  Eqs.  (lU)  through  (l6).  A  more  complete  discussion  is  given 
in  Appendix  A.  When  9^0,  Eq.  (15)  is  no  longer  the  Mathieu  equation, 
and  while  Eq.  (15)  admits  a  series  solution  of  a  type  similar  to  the 
Mathieu  functions,  a  search  of  the  literature  has  failed  to  locate 
any  treatment  in  detail  of  the  solution  or  characteristic  numbers 
associated  with  Eq.  (15).  Therefore,  as  shown  in  Appendix  A,  all  we 
can  say  about  Eq.  (15)  without  going  into  the  details  of  a  series 
solution  is  that  the  solution  w  -  A^w^  +  RjW^  is  of  the  form  where 
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K,(t) 

wi  "  e  1  (26) 

where  the  characteristic  numbers  ^  are  related  to  the  roots  of  the 
quadratic  characteristic  equation  (Eq.  (A-IO))  and,  as  functions  of  the 
free  parameters  0  and  0,  admit  a  complete  spectrum  of  real,  imaginary, 
and  complex  values.  The  functions  ^(t)  are  periodic,  of  period  2n, 
and  are  hence  bounded.  Therefore,  the  solutions  v.^  are  bounded  for 
all  t  only  when  both  characteristic  numbers  are  imaginary.  Let  us 
write  these  solutions  in  the  expression  for  2  given  by  Eq.  (l4): 


f  H_(t)  1  icp  -3t/2 

-  •<  $i(r )  +  Age  §2(t)  >•  e  e  (27) 


It  Is  immediately  obvious  from  Eq.  (27)  that  |  z|  will  be  bounded  for 
all  time  t  if  and  only  if 


I  Re  (u±)|  <  5/2  (28) 

mr 

That  is,  if  ui  is  complex,  the  aperiodic  portion  of  e  will  cause 
|z|  to  become  infinite  unless  the  damping  factor  e  predominates. 

Let  us  examine  this  criterion  in  comparison  with  the  solution  of 
the  previous  section  where  3-0.  There,  the  instability  of  w1  or  w^ 
implied  immediately'  that  |z|  would  be  unbounded  as  t  -»  because  the 
damping  term  e~^T'2  was  lacking  in  the  solution.  Here,  we  have  intro¬ 
duced  a  term  which  would  apparently  increase  the  stability  of  the 
system.  It  is  intuitively  obvious,  but  difficult  to  prove  without  a 
series  solution,  that  this  is  indeed  the  case.  For,  since  n ^  will 
depend  on  both  3  anu  Q,  the  regions  of  stability  of  Eq.  (15)  may  be 
less  extensive  than  those  of  the  Mathieu  equation,  and  Eq.  (27)  might 
have  regions  of  stability  which  are  less  extensive  than  those  of  Eq. 
(17)  and  (18). 
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Therefore,  the  only  definite  conclusion  that  can  be  made  about 
Eq.  (27)  is  that  stable  solutions  exist.  If  the  characteristic  num¬ 
bers  are  such  that  ^  and  are  complex  conjugate,  and  if  Re(u^) 

*  0/2  then  jz|  is  a  periodic  function  of  t.  Otherwise,  when  |Re(u^)| 

<  3/2,  |  z|  will  approach  zero  as  t  -*  »,  (This  is  true  when  the  char¬ 
acteristic  numbers  ^  are  imaginary  as  well. ) 

The  preceding  discussion  has  been  for  arbitrary  0  and  Ci.  How¬ 
ever,  the  averaging  process  described  earlier  demands  that  certain 
restrictions  on  3  and  Q,  given  by  Eqs.  (10)  and  (13),  be  met.  Within 
this  structure,  it  is  possible  to  average  the  motion  over  times  long 
compared  to  the  cyclotron  frequency  if  Q  »  1  and  to  obtain  the  motion 
of  the  guiding  center;  hence  one  case  of  interest  will  be  that  in 
which  Cl/ 3  <  1,  0  »  1.  On  the  other  hand,  if  one  averages  over  times 
short  compared  to  the  cyclotron  period,  then  Cl  is  not  restricted,  and 
Eqs.  (10)  and  (13)  combine  to  give  Cl/ 3  «  1,  Q  ^  l, 

AN  APPROXIMATE  SOLUTION 

The  fact  that  3  »  1  leads  one  to  consider  an  approximate  solu¬ 
tion  cf  Eq.  (15)*  Rewrite  Eq.  (15)  as 


V”  +e2  (7.  -Si’Ll  tl)2  w  .  0  (29) 

where  7  ■  0/3  is  the  ratio  of  the  cyclotron  frequency  to  the  collision 
frequency.  From  Eq.  (lO),  7  ~  1  when  the  Langevin  time  average  is 
taken  over  a  cyclotron  period;  the  second  case,  in  which  the  time 
average  is  much  less  than  a  cyclotron  period,  is  given  by  Eq.  (13), 
and  7  «  1. 

The  method  of  solution  when  $  »  1  (sometimes  called  the  WKB 
method)  is  discussed  fully  In  Appendix  B.  The  results  are  reproduced 
here  for  convenience.  The  general  solution  of  Eq.  (29)  when  terms  of 
the  order  l/32  are  neglected  is 
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w 


-/l  -  i y  sin  t 


Icp-St/2  -Icp+Pt/2- 


+  Age 


(30) 


where  cp  «  fi  sin^T/2. 

It  is  obvious  from  the  undamped  exponential  in  the  second  term 
that  w  is  not  bounded  as  t  -*  «.  We  should  note  also  that  v(t)  has 
the  general  appearance  predicted  by  Eg.  (26)  of  an  exponential  times 
a  periodic  function  of  t. 

The  fact  that  v(t)  is  unbounded  does  not  preclude  a  bounded 
solution  of  Eg.  (27);  a  comparison  of  the  exponential  terms  in  Eqs. 
(27)  and  (30)  shows  that  the  undamped  exponent  will  vanish  when  the 
two  equations  are  combined,  viz. 


z 


(1 


2io 


(31) 


where  the  initial  conditions  on  z  have  been  used  to  determine  the  con¬ 
stants  of  integration  and  A^.  After  considerable  algebra,  z  can 
be  written  in  camparitively  simple  form: 


z 


6(1  ♦  r2  sin2  t)1'1*  ^ 


-  F 

|_|  ten  1(y  sin  t)  +  i  ten"1  jrj 


(32) 


where 


F1  *  V  +  0  -  e  cos  2c  -  (C^  -  sin  (33) 


-Pt 


L(0i  -  I> 


cos  &?  +  V  sin 


ap] 


(34) 


-2U- 


remembering  that  cp  ■  ft  sin^T/2,  and  y  ■  n/P. 
The  components  of  motion  are  thus 


I2  - 


*i +  i 


P2(l  +  72  sin  t) 


7F* 


(35) 


and 


.  ,  f(F?  -  F?)  7  sin  x  +  2F.Fr 

arg  z  ■  —  ten  J - - 


1(1^  -  F^)  -  2F 


■n  x  *  2F1Fg  1 
'^F^  sin  t  J 


(36) 


It  is  apparent  from  Eqs.  (35)  and  (36)  that  the  orbit  will  be 
bounded  for  all  time  t  >  0,  and  if 


+  ^ 

$ — ~  <  R0  (37) 

where  R  is  the  dimensionless  radius  of  the  container,  that  the  motion 
o 

will  be  confined  to  regions  inside  the  tube  for  t  >  0.  This  is  true 
regardless  of  the  stability  conditions  on  Eq.  (29)>  since  the  damping 
term  e’PT  removes  the  instability  present  in  Eq.  ( 30). 

Although  Eqs.  (35)  and  (36)  are  complicated,  the  motion  in  reality 
is  quite  simple.  Since  P  »  1,  the  damping  terms  in  F1  and  Fg  will 
vanish  for  relatively  small  times,  and  if  we  further  assume  that 
V  *  ■  0  suad  7*1,  then 

R2  „  1  ±  72AP2  w  —-±-— 

7l  +  y2  sin2  t  +  72  sin2  t 


and 


(38) 
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when  7/P  *  O/P  £  1,  This  shows  that  the  motion  is  a  simple  oscilla¬ 
tion  in  which  the  electron  position  varies  in  a  periodic  manner,  reach¬ 
ing  its  maximum  extent  when  the  electric  field  is  zero.  Highly  damped 
motion  corresponds  to  ohmic  resistivity,  as  would  be  expected.  The 
maximum  extent  of  the  motion  is  at  a  radial  position  of  R  »  0.84,  and 
the  maximum  angular  displacement  is  7  radians. 

One  would  suspect  from  Eqs.  (38)  and  (39)  that  the  electron  is 
experiencing  very  little  acceleration  in  its  trajectory.  In  fact,  if 
one  were  to  solve  the  original  differential  equation,  Eq.  (7),  with 
the  acceleration  term  z’!  »  0,  the  result  would  be  identical  with  the 
approximate  versions  of  Eqs.  (37)  and  (38),  in  which  terms  of  the 
order  7/P  are  neglected.  This  result  is  not  unexpected  from  the 
Lange vln  equation,  because  the  drag  terra  there  is  proportional  to  the 
velocity.  When  the  collisions  are  predominant,  the  particle  has  no 
chance  to  accelerate  before  the  drag  decelerates  it  again. 

NUMERICAL  EXAMPLES 

A  few  trajectories  have  been  numerically  computed  to  illustrate 
the  motion  of  an  average  electron  when  collisions  predominate.  The 
values  of  0  chosen  for  these  calculations  are  the  same  as  those  in 
the  collisionless  motion.  The  parameter  7  (the  ratio  of  cyclotron 
frequency  to  collision  frequency)  was  varied  from  7  ■  1  to  7  ■  l/4. 
These  results  are  presented  in  Figs.  7  and  8. 

For  purposes  of  comparison,  one  value  of  0  was  selected  so  that 
when  3*0  the  trajectory  became  identical  with  that  computed  for  the 
collisionless  case.  Figure  7,  for  example,  shows  the  motion  for  0  * 
17*242  and  various  values  of  3;  the  collisionless  motion  for  this  case 
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Average  motion  of  electrons  for  Cl  ■  17. 242 
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Fig. 8 — Average  motion  of  electrons  for  il=  2.316 
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is  given  in  Fig.  3.  TCie  unstable  collisionless  motion  for  0  ®  2.316 
becomes  stable  vhen  0  »  1.  An  example  of  this  is  given  in  Fig.  8. 

ttie  motion  given  in  Figs.  7  and  8  agrees  well  with  the  approxi¬ 
mate  solution  given  by  Eqs.  (38)  and  (39)  when  t  is  large  enou^i  so 
that  e~^T  -*  0.  The  complete  approximate  solution  given  by  Eqs.  (35) 
and  (36)  should,  in  fact,  more  than  adequately  describe  the  motion  as 
long  as  3  is  larger  than  the  limits  prescribed  in  Appendix  B. 

In  conclusion,  one  can  say  that  the  average  motion  of  an  electron 
in  a  weakly  ionized  gas  is  bounded  in  the  sense  that  the  position  of 
the  electron  oscillates  90  deg  out  of  phase  with  the  applied  field 
and  will  not  intersect  the  container  wall  if  Eq.  (37)  is  satisfied. 
However,  since  diffusion  is  neglected,  some  electrons  will  always 
leave  the  high-density  region  for  the  wall,  and  thus  the  mathematical 
existence  of  bounded  solutions  does  not  imply  confinement. 
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Appendix  A 

ON  SECOND- ORDER  LINEAR  DIFFERENTIAL  EQUATIONS 
WITH  PERIODIC  COEFFICIENTS 

The  basic  equation  of  this  study,  Eq.  (15),  is  of  the  form 


—  +  p(t)v  -  0  (A-l) 
dr 

where  p(r)  is  a  continuous,  periodic,  complex  function  of  the  real 
variable  r  over  the  range  -  «  <  r  <  ».  The  function  p(t)  contains 
two  general  parameters,  Q  and  0. 

The  most  familiar  form  taken  by  Eq.  (A-l)  is  that  for  the  case 
in  which  p(t)  Is  a  real,  even,  periodic  function  P(t);  Eq.  (A-l)  is 
then  generally  called  the  Hill  equation: 


&  +  (P(t  )  +  X)w  -  0  (A-2) 

dT 

where  X  is  some  areal  parameter.  Heare  t  can  be  either  areal  or  complex. 
A  special  case  of  the  Hill  equation  is  the  Mathieu  equation 


♦  (a  -  2q  cos  2r)  w  -  0  (A-3) 

d tZ 

where  a  and  q  a  are  again  paarameters  in  the  equation. 

The  theory  of  the  genearal  solution  and  its  behavior  for  the  case 

in  which  p(t)  is  aneal,  such  as  in  Eq.  (A-2)  or  (A-3),  is  discussed  at 

some  length  in  Ince^®^  and  Mela ch lan.  The  fact  that  the  paresent 

equation  has  the  same  type  of  genearal  solution  was  not  known  to  the 

(19) 

author  until  it  was  found  in  Coddington  and  Levinson.  Because 
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the  nature  or  stability  of  the  general  solution  of  Eq.  (A-l)  is  essen¬ 
tial  to  the  discussion  in  the  text,  a  short  discussion  of  its  existence 
and  properties  is  Included  here* 

It  can  be  shown  (theorems  2*3  -  5*1  of  Ref-  18)  that  Eq.  (A-l) 
admits  a  solution  of  the  form 

M*i T  kUT 

v(t)  -  A^e  X  #1(t )  +  A^e  d  $2(t)  (A-U) 

where  *(t )  is  periodic,  of  the  same  period  of  p(t),  and  where  the  ex¬ 
ponent  u  is  determined  by  the  roots  of  the  characteristic  equation  and 
is  a  function  of  the  parameters  0  and  3.  An  outline  of  the  proof  is 
as  follows:  The  set  of  all  solutions  of  Eq.  (A-l)  forms  a  two-dimen¬ 
sional  vector  space  over  the  complex  field.  Two  of  these  solutions 
are  linearly  Independent  and  hence  form  a  basis  for  the  space*  These 
are  called  the  fundamental  set.  Any  linear  combination  of  these  is 
also  a  solution.  If  v^(t)  and  v2(t)  are  the  two  solutions,  then 

v(t)  ■  a3w1(t)  +  \v2(t)  (A-5) 

la  a  solution.  Now  p(t  +  2*1 )  ■  p(r),  since  2"  is  the  period  of  p(t) 
in  this  case.  Therefore  v^(t  +  2 n)  and  Wg(r  +  2t)  satisfy  Eq.  (A-l) 
and  can  also  be  written  in  terms  of  the  fundamental  set  as 


w^t  ♦  2rr)  -  OjjV^t)  +<>jL2w2(T) 
vg(T  ♦  2»r)  -  c^w^t)  +  o^v^t) 


(A -6) 


Writing  Eq.  (A-5)  at  r  ♦  2n,  and  coariblnlng  this  with  Eq.  (A-6),  one 
has,  for  the  general  solution  at  t  + 


w(t  +  2")  ■  (AjOt^  +\a^1W1(T)  +  (Aj0^  +\°fe2^T^ 
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i.e,,  w(t  +  2n)  is  a  linear  combination  of  the  fundamental  solutions 
v^(t)  and  w^t).  We  want  to  determine  a  general  solution  w(t)  such 
that 

v(t  +  2n)  ■  kw(T)  (4-8) 

This  Is  possible  if  k  is  related  to  the  coefficients  in  Eq.  (A-7)  in 
the  following  way: 

^3  *  Afn  +  V21 

(A-9) 

*  A3«L2  +  V22 

which  have  a  solution  if  and  only  if 

°L1-  k  °fel 

=■  0  (A  10) 

“12 

The  deteralnantal  Eq.  (A-10)  is  called  the  characteristic  equation. 
and  k  the  characteristic  roots  of  this  equation.  If  we  let  k  ■  e^1, 
then  u  is  the  characteristic  exponent,  which  is  determined  to  within 
modulo  2n  for  the  imaginary  part.  The  real  part  of  n  is  determined 
uniquely. 

Let  us  write  the  general  solution  as 

w(t)  -  eMT*(T)  (A-ll) 

where  *(t)  Is  unspecified  for  the  moment.  Writing  Eq.  (A-ll)  at  t  +  2n, 

we  have 

#(t  +  2n)  ■  e"^T_2rT^  w(t  +  2")  (A-12) 
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and,  using  Eq.  (A-8)  and  the  definition  of  k 

$(t  +  2n)  ■  e~UTw(T)  =  ®(t)  (A-13) 

Therefore  $(t)  is  a  periodic  function  of  t  with  the  same  period  (2tt) 
as  the  coefficient  p(t). 

To  conclude,  we  have  shown  that  Eq.  (A-l)  has  the  solutions 

*  e  1  *1(t)  (i  ■  1,2)  (A-lU) 

where  the  two  characteristic  exponents  are  given  by  the  solution 
of  Eq.  (A-10).  The  coefficients  a ^  in  Eq.  (A-10)  are  found  from  the 
initial  data  and  the  solution  evaluated  at  t  -  2n  and  are  hence  func¬ 
tions  of  the  parameters  0  and  3.  The  general  solution  of  Eq.  (A-l) 
is  therefore  given  by  Eq.  (A-4). 

THE  NATURE  OF  THE  SOLUTION 

First,  we  can  write  the  coefficients  of  the  characteristic  equa¬ 
tion  in  terms  of  the  initial  conditions  and  the  solutions.  Let  us 
pick 


w1(0)  -  1  v2(0)  -  0 

v£(0)  -  0  Wg(o)  -  1 


(A-15) 


and  let  the  constants  A^,  in  Eq.  (A-4)  take  care  of  the  initial 
conditions  on  Eq.  (lj)  of  the  text.  Then  Eq.  (A -6),  differentiated 
and  evaluated  at  r  *  0,  yields 

Vi(2»T)  -  Oj^ 


v'(2n)  -  0^2 


(A-16) 
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Expanding  the  determinant  given  toy  Eq.  (A-10),  we  have 

k2  -  ^(2")  +  v^(2tt)^  k  +  1  -  0  (A-17) 

where  use  has  been  made  of  the  fact  that  the  Wronskian  of  the  funda¬ 
mental  set  v^w^  -  W2W1  "  Equation  (A-17)  shows  that  the  values  of 
k  depend  on  the  values  of  the  solution  at  the  point  t  *  2rr.  When  k 
Is  such  that  y  is  purely  imaginary,  the  solution  remains  finite  as 
r  -*  ».  On  the  other  hand,  if  y  is  real  or  complex,  then  e4T  becomes 
infinite  for  either  t  »  or  t  In  this  case  the  solution  is 

unstable. 

While  the  theory  Just  presented  is  very  simple,  it  is  not  easy 

(19) 

to  determine  the  characteristic  exponents.  Ince  discusses  a 
method  for  doing  so  which  involves  repeated  integration  of  p(t)  with 
the  terms  of  a  power  series  expansion  of  the  fundamental  set  of  solu¬ 
tions.  These  solutions  could  be  obtained  by  a  series  expansion.  Be¬ 
cause  of  the  nature  of  p(t),  however,  the  undertaking  of  such  a  task 
is  not  warranted.  Even  with  the  simple  Mathieu  equation,  a  whole 
book^^  is  devoted  to  the  various  series  solutions,  stability  dia¬ 
grams,  and  numerical  techniques. 

THE  MATHIEU  EQUATION 

The  Mathieu  equation  emerges  from  Eq.  (A-l)  when  0*0: 

+  (jj—  sin^  t)  w  *  (l  -  cos  2t)  w  *  0  (A-l8) 

dt  dT 

Notice  that  now  p(t  )  is  an  even  function  of  t  and  permits  us  to  write 
for  a  general  solution 


w(t)  *  Ae^T*(T)  +  Be~MT$(-T) 


(A-19) 
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Let  ?(±t)  »  Z  c  e±iJDT  and  let  n  ■  K  +  la.  Then 

-oo  m 


v(t) 


*  r 

V  -5 


c  el(crhn)T 

IQ 


+  Age 


-Kt 


5  c  e'l(<rtn)T 


(A-20) 


This  rather  simple  form  for  the  series  solution  is  due  to  the  even¬ 
ness  and  periodicity  of  p(r).  In  the  general  case  (Eq.  (A-l)),  a 
series  solution  would  be  much  more  complicated. 

If  li  *  ia,  when  a  is  nonintegral,  Eq.  (A-20)  is  stable,  (if  a 
is  integral,  the  solution  is  called  neutral. )  If  a  is  rational,  say 
a  m  p/s,  then  the  period  of  w  is  2ns.  If  a  is  irrational,  w  is  aperi¬ 
odic;  i.e.,  oscillatory  and  bounded,  but  nonrepeating. 

Ihe  characteristic  roots  of  Eq.  (A-10)  have  been  found  for  this 
set  of  fundamental  solutions,  which  are  Mathieu  functions.  Figure  9 
reproduces  in  part  the  regions  of  stability  as  a  function  of  the  para¬ 
meter  0,  along  with  some  points  indicating  those  values  of  ft  which 
yield  values  of  a  giving  period  of  4n.  The  unstable  solution  illus¬ 
trated  in  Fig.  6  is  also  indicated  in  Fig.  9*  The  legend  "even"  and 
"odd”  in  Fig.  10  (Appendix  B)  indicates  whether  the  index  in  Eq.  (A-20) 
is  even  or  odd. 

The  numerical  determination  of  ft  for  a  m  l/2  so  that  ft  lies  in 
one  of  the  stable  regions  is  treated  extensively  in  Ref.  15.  The 
author  found  that  the  method  of  continued  fractions  (8  5. 11  of  Ref. 

15)  was  the  simplest  to  perform. 

The  lines  dividing  regions  of  stable  and  unstable  solutions  cor¬ 
respond  to  Mathieu  functions  of  integral  order.  Here  a  ■  0  or  a  ■  1, 

(18I 

and  it  is  not  possible  to  have  two  periodic  solutions.  For  that 

reason,  the  examples  given  in  the  text  are  for  fractional  values  of  a. 
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Type  of  solution 


Fig.  9  —  Stability  chart  for  solutions  of  Mathieu 
equation,  showing  values  of  A  for  which 
numerical  solutions  were  obtained 
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Appendix  B 


APPROXIMATE  SOLUTION  OF  Eg.  (15) 


Write  Eq.  (15)  as 


&  ♦ (i  as.  x  1  i)\ .  o 

djd  * 


(B-l) 


vhere  7  ■  z  <  and  0  »  1.  Let  x  ■  v'/Bv  and  p  ■  (7  sin  t  +  i)/2; 


then 


1  dx  .  2  .  2 
0  dx  x  9 


(B-2) 


(This  is  the  Riccati  equation. )  Assume  that 


.  1  ,  1  . 

Xo  3X1  g2  X2  #  # 


(b-3) 


Inserting  this  expression  for  x  into  Eq.  (B-2)  and  equating  the  co¬ 
efficients  of  powers  of  3  to  zero  gives 


x  .  ±  10  .  Sl  ,  L.  (2PP11  -  S'x2  + 

X  *  *  1P  20P  *  g2  (  ^3  '  •" 


(B-t) 


so  that 


.Am  ;  s  4  app"-|'.'!  a.  ♦ ... 

V  p  8p3 


(B-5) 


Neglecting  the  integral  term  for  the  time  being,  we  note  that  v(t) 
as  given  by  Eq.  (B-5 )  is  unstable,  since  the  term  e®T^2  -•  «  as  t  -»  ». 
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However,  the  motion  itself  is  stable,  which  is  seen  if  ve  substitute 
w(t)  into  Eq.  (14),  yielding 


z 


1  +  ±7  sin  t 
1+7^  sin^  t 


|(V  +  iCx  -  |*)(l 


e2icp-Pt 


(B-6) 


This  is  Eq.  (31)  of  the  text. 

The  neglect  of  the  integral  term  introduces  an  error  which  is 
best  seen  by  examining  Eq.  (B-4).  Since  p  *  (7  sin  t  +  i)/2,  we  can 
write  the  term  containing  p ,  p  1 ,  and  p  '  *  in  terms  of  p .  Letting  the 
last  term  in  Eq.  (B-4)  be  given  by  f(p),  then 


w 


[iRe(p  -  f(p))]  -  [Im(p  -  f(p))]ds 


(B-7) 


where 


Re(p  -  f(p)) 


y[y4  -  (9  +  372)y2  +  (9r2  +  2)] 

2  ’  402(y2  +  l)3 


(B-8) 


Im(p  -  f(p)) 


1  +  5y4  -  (7  *  972)y2  ♦  372 

2  402(y2  +  l)3 


(B-9) 


and  where  y  «  2p  -  i  «  7  sin  t  is  veal*  The  contribution  of  the  second 
terms  to  the  integrand  is  largest  (for  all  time)  when  7  ~  1*  When 
7  «  1,  then  y  «  1,  and  the  integrand  becomes 

1  •  <b-i°> 


in  which  the  dominant  term  Is  the  real  part. 
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computing  values  of  Eqs.  (B-8)  and  (B-9)  for  7  ■  1,  it  is 
possible  to  estimate  the  lower  limit  on  6  below  which  the  approximate 

p 

solution  is  invalid.  Figure  10  shows  the  two  equations  for  0  *  1, 

10,  and  «  over  a  quarter-period  of  the  function  f(p).  It  is  obvious 
that  the  error  induced  in  the  periodic  portion  of  the  solution  is 

o 

small  when  0  >  10.  The  error  in  the  real  part  of  the  integrand  will 

p 

average  out  over  a  period  of  f(p)  for  0  >  10.  Hence  0  >  3  would  ap¬ 
pear  to  satisfy  the  approximation  in  Eq.  (B-6). 

It  is  also  interesting  to  note  that  the  approximate  solution  given 
by  Eq.  (B-j),  when  substituted  back  into  the  differential  equation 
Eq.  (B-l),  satisfies  Eq.  (B-l)  whenever  0^  >  10  (or  f(p)  **  0).  This 
is  true  whether  the  integral  term  in  Eq.  (B~5)  is  retained  or  not. 


Fig.  10 — Estimate  of  the  error  in 
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